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1.  Introduction 


Reactive  materials  can  be  loosely  categorized  as  composites  of  inert  solid  materials  which,  when 
subjected  to  a  violent  mechanical  stimulus  such  as  an  impact,  react  exothermally  with  a  rapid 
release  of  energy.  This  reaction,  while  aptly  described  as  “explosive,”  differs  from  a  true 
detonation  or  deflagration  in  that  it  requires  a  mechanical  stimulus  to  not  only  initiate  the 
reaction  but  also  to  sustain  it.  Such  materials  can  also  be  fairly  robust  mechanically  and  can 
serve  as  substantial  structural  components.  Because  of  these  properties,  reactive  materials  have  a 
number  of  potential  ordnance  applications.  Various  compositions  have  been  investigated  to 
tailor  properties  of  reactivity,  strength,  and  density,  to  suit  particular  needs. 

This  report  addresses  a  specific  aluminum/polytetrafluoroethylene  (Al/PTFE)  formulation  that 
serves  as  an  important  benchmark  for  current  reactive  material  development.  In  addition  to 
material  development,  a  good  deal  of  work  is  underway  to  develop  physics-based  modeling 
capabilities  for  reactive  materials.  To  aid  in  this  latter  effort,  this  report  provides  basic 
constitutive  models  for  the  inert  behavior  of  this  material.  Compression  tests  have  been 
performed  over  a  range  of  strain  rates  and  temperatures  relevant  to  the  conditions  present  during 
low-speed  impact.  This  data,  presented  in  the  following  sections,  is  used  to  generate  parameters 
for  both  the  Johnson-Cook  (JC)  and  Modified  Johnson-Cook  (MJC)  constitutive  models. 
Although  not  ideally  suited  to  represent  this  material,  these  parameters  are  given  primarily  due  to 
the  wide-spread  use  of  these  models  and  their  availability  in  the  current  suite  of  hydrocodes.  An 
additional,  and  more  appropriate,  fit  is  given  for  the  Zerilli-Armstrong  (ZA)  model  for  polymers. 

The  remainder  of  this  document  is  organized  as  follows:  a  basic  material  description  is  given  in 
section  2;  section  3  describes  the  testing  techniques  used  in  the  experiments,  and  the  data  is 
presented  in  section  4;  finally,  the  model  parameters  and  descriptions  of  the  fitting  processes  are 
given  in  section  5.  Unless  otherwise  noted,  the  use  of  the  words  “stress”  and  “strain”  denote 
positive  values  in  compression. 


2.  Material 


The  samples  tested  during  this  program  were  supplied  by  General  Sciences,  Inc.  (GSI).  They 
were  made  from  a  material  designated  by  GSI  as  GSI-00I7.  It  is  a  pressed  and  sintered  mixture 
of  aluminum  and  PTFE  powders,  26.5%  and  73.5%  by  weight,  respectively.  The  initial  powder 
sizes  are  44  and  3 1  |im,  respectively.  The  density  of  the  compacted  material,  as  measured  by  a 
buoyancy  method  based  on  Archimedes  principle,  is  2.29  g/cm  .  The  longitudinal  wave  speed, 
measured  by  ultrasound  at  10  MHz,  is  1570  m/s. 
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3.  Experimental  Procedure 


3.1  Low-Rate  Tests 

Low-rate  tests  were  performed  with  an  Instron,  Ine.,  model  1331  servo-hydraulie  load  frame. 

The  load  applied  to  the  speeimen  was  measured  with  a  load  eell,  and  the  speeimen  deformation 
was  measured  using  a  linear  variable  differential  transformer  (LVDT)  measurement  of  the  eross- 
head  displacement  and  includes  a  correction  for  machine  compliance.  The  specimens  were 
cylindrical,  nominally  6.35  mm  in  both  diameter  and  length.  Contact  surfaces  were  lubricated 
with  a  heat-stable  silicone  lubricant.  All  of  the  low-rate  tests  were  performed  at  room 
temperature  (22  °C). 

3.2  High-Rate  Tests 

A  6.35-mm-diameter  7075-T6  aluminum  Split  Hopkinson  Pressure  Bar  (SHPB)  was  used  for  the 
high-rate  tests.  A  series  of  six  tests  were  performed  at  rates  from  600  to  8000/s,  all  initially  at 
room  temperature  (22  °C).  The  specimens  were  cylindrical,  3.18  mm  in  diameter  and  length, 
and  contact  surfaces  were  lubricated  with  the  same  silicone  lubricant  used  in  the  low-rate  tests. 
All  data  analysis  accounts  for  bar-wave  dispersion  using  a  method  based  on  Follansbee  and 
Franz  (7)  and  Gong  et  al.  (2). 

As  is  typical  with  the  SHPB,  axial  deformation  of  the  specimen  is  determined  by  the  bar  analysis 
(see,  for  example,  Follansbee  [5]).  However,  the  radial  deformation  is  unknown,  and  is  usually 
determined  under  the  assumption  that  the  deformation  (both  elastic  and  inelastic)  in  the  specimen 
occurs  at  constant  volume.  To  verify  this  assumption,  a  laser  extensometer  (also  known  as  a 
laser  occlusive  radius  detector,  or  LORD  [4])  was  used  to  measure  the  radial  deformation  of  the 
specimens  during  dynamic  compression.  The  frequency  response  of  this  device  is  very  high 
(-'10  MHz),  and  in  our  experience  it  is  able  to  resolve  deformations  as  small  as  0.01  mm.  The 
results  of  these  experiments  are  discussed  in  section  4. 

3.3  Temperature  Tests 

A  final  set  of  experiments  was  performed  at  elevated  temperatures  to  quantify  the  thermal 
softening  behavior  of  the  material.  These  were  performed  at  a  consistent  strain  rate  of  4000/s 
using  the  SHPB.  Heating  was  accomplished  by  circulating  heated  air  into  a  chamber  which 
enclosed  the  specimen  and  the  adjacent  ~65  mm  sections  of  the  bars.  Ideally,  specimen 
temperature  would  have  been  monitored  directly  with  a  thermocouple  glued  to  each  specimen. 
However,  this  proved  impractical  because  of  the  small  sample  size  and  also  because  of  the 
difficulty  in  adhering  gages  to  the  specimen.*  Instead,  specimen  temperature  was  measured  with 
a  thermocouple  probe  placed  within  1  cm  of  the  specimen,  i.e.,  the  probe  measures  ambient  air 


PTFE,  more  commonly  known  as  Teflon,  is  widely  known  for  its  “non-stick”  properties. 
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temperature  and  not  the  specimen  temperature  directly.  The  temperature  in  the  chamber  was 
allowed  to  equilibrate  over  a  20  min  period  prior  to  each  test  to  ensure  that  the  specimen  and 
relevant  sections  of  the  bars  were  allowed  to  reach  the  ambient  temperature.  The  temperature 
gradient  in  the  bars  is  believed  to  have  negligible  effects  on  the  bar  wave  propagation  and  the 
strain  gage  measurements.  Temperature  measurements  made  in  this  way  are  estimated  to  be 
accurate  to  within  ±2  °C. 


4.  Results 


Figure  1  shows  the  engineering  stress-strain  curves  obtained  from  the  room  temperature  tests. 
Because  of  the  brief  duration  of  the  deformation,  the  high  rate  curves  (600/s  and  beyond)  are 
considered  adiabatic.  In  contrast,  the  lowest  rate  curves,  0.001  and  0.01/s,  are  considered 
isothermal.  The  intermediate  rate,  at  0.1/s,  is  probably  somewhere  between  the  two  limits. 


Figure  1.  Engineering  stress  vs.  engineering  strain  for  Al/PTFE  (compression).  All  specimens  were 
initially  at  room  temperature  (22  °C). 

Figure  2  shows  the  tensile  radial  engineering  strain  as  measured  by  the  laser  extensometer  as  a 
function  of  compressive  axial  engineering  strain  for  the  six  SFIPB  experiments  shown  in 
figure  1 .  These  are  plotted  as  thin  black  curves.  The  heavy  red  curve  is  the  relationship  that 
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Figure  2.  Lateral  strain  as  measured  by  a  laser  extensometer  in  the  SFIPB  tests  plotted  in  figure  1  (thin 
black  curves).  The  heavy  red  curve  shows  that  which  would  be  expected  from  incompressible 
behavior. 

would  result  from  incompressible  deformation.  The  agreement  is  good,  and  what  deviation  can 
be  measured  could  easily  be  explained  by  barreling.  Therefore,  the  remainder  of  this  report  will 
speak  in  terms  of  true  stress  and  strain  as  calculated  from  the  engineering  values  under  the 
assumption  of  incompressibility.  These  terms  and  the  assumption  of  incompressibility  are 
reviewed  in  the  appendix. 

Figures  3  and  4  show  the  true  stress-strain  curves  for  the  tests  of  figure  1 .  The  data  are  plotted 
with  two  scales  to  show  more  detail  in  the  low-strain  range.  Figure  5  shows  the  rate  of  true 
strain  for  the  SFIPB  experiments.  Note  that  in  all  cases  the  strain  rate  is  reasonably  constant. 

The  low-rate  tests  are  controlled  in  the  servo-hydraulic  load  frame  to  give  constant  rates  of  true 
strain.  Figure  6  shows  stress  at  7%  strain  as  a  function  of  strain  rate,  and  shows  the  trend  of  rate¬ 
hardening  of  this  material. 

Figure  7  shows  the  stress-strain  curves  from  the  experiments  at  4000/s,  which  were  conducted  at 
initial  temperatures  of  22,  43,  56,  and  78  °C.  A  substantial  softening  is  observed  over  this  range 
of  temperatures.  Stresses  at  a  strain  of  10%  are  plotted  with  temperature  in  figure  8.  From  this 
graph,  it  is  apparent  that  the  softening  is  approximately  linear  with  stress  over  the  range  of 
22-56  °C.  As  seen  in  figure  7,  there  is  a  significant  change  in  the  form  of  the  stress-strain  curve 
for  the  78  °C  test.  It  is  suggested  that  this  may  be  due  to  a  phase  change  in  the  PTFE  that  occurs 
at  some  temperature  between  56  and  78  °C,  although  this  is  merely  speculative. 
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Figure  3.  True  stress-strain  eurves  for  AFPTFE.  All  samples  were  initially  at  room  temperature  (22  °C). 


Figure  4.  The  data  of  figure  3  repeated  on  a  finer  seale  to  illustrate  the  low-strain  behavior. 
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Figure  5.  The  variation  of  true  strain  rate  for  the  SFiPB  tests. 


strain-rate  (1/s) 


Figure  6.  Strain-rate  hardening  observed  in  the  experiments  compared  to  the  MJC,  JC,  and  ZA  fits. 
Stress  is  measured  at  7%  total  true  strain,  all  at  room  temperature  (295  K). 
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true  strain 


Figure  7.  Stress-strain  curves  obtained  at  a  rate  of  4000/s  over  a  range  of  initial  temperatures. 


temperature  (K) 


Figure  8.  Three  fits  to  the  high-rate  temperature  data  (10%  total  strain,  4000/s).  The  JC  equation  is 
fit  over  two  temperature  ranges:  the  lower  temperature  range  of  295  K  <  T  <  329  K  and 
the  full  temperature  range.  A  fit  is  also  shown  for  the  ZA  model  (full  range). 
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5.  Constitutive  Modeling 


5,1  Johnson-Cook  and  Modified  Johnson  Cook 

The  data  of  the  preeeding  seetions  have  been  fitted  to  the  JC  and  MJC  eonstitutive  equations. 
These  models  were  ehosen  primarily  beeause  of  their  widespread  use,  and  not  neeessarily  beeause 
they  are  partieularly  well-suited  to  this  material.  Details  of  the  models  are  diseussed  in  Johnson 
and  Cook  (5)  and  Holmquist  and  Johnson  (6),  and  are  reviewed  briefly  in  the  following. 

The  JC  and  MJC  equations  are  written  as 


and 


C)=(A  +  B£p’')  1  +  Cln  —  (l-T*”) 


(1) 

(2) 


respeetively,  where  T*  is  the  homologous  temperature, 

T-T 

T*  = - (3) 

The  term  £p  denotes  plastie  strain,  i.e.,  the  total  strain  less  the  elastie  strain,  assumed  in  the 
following  fits  to  be  3.5%.  A,  B,  n,  C,  A,,  and  m  are  all  material  eonstants.  Tm,  whieh  is  generally 
taken  as  the  melting  temperature,  is  here  treated  as  a  free-parameter  to  improve  the  quality  of  the 
fits.  Tr  is  a  referenee  temperature,  295  K,  the  temperature  of  the  lowest  temperature  tests.  Note 
that  the  form  of  the  thermal  softening  term  is  in  general  not  defined  for  T  <  Tr,  so  extrapolation 
below  Tr  is  not  possible.  The  referenee  strain-rate.  Eg ,  is  taken  as  1/s. 

The  only  differenee  in  the  two  equations  is  the  form  of  the  strain-rate  hardening;  the  original 
version  uses  a  logarithmie  form  and  the  modified  version  uses  a  power-law  form.  The  form  of 
the  equations,  a  =  f(£)-g(£)-h(T),is  very  eonvenient  for  fitting  the  available  data.  By 
assuming  £  ~  0  and  T  =  295  K  in  equation  1 ,  estimates  of  A  and  C  ean  be  determined  from  a 
simple  least  squares  fit  to  the  rate  data  of  figure  6.  This  ean  then  be  repeated  with  equation  2, 
giving  estimates  for  A  and  A  for  the  MJC  equation. 

Next,  the  temperature  data  of  figure  8  ean  be  fit  to  eaeh  equation  under  the  assumptions  that 
£  ~  0  and  £  =  4000/s  to  give  estimates  of  A,  m,  and  Tm.  Two  fits  are  made,  one  over  the  lower 
three  temperatures,  a  range  of  T  G  [295, 329]  K,  and  the  other  the  full  range  T  G  [295,351]K.  The 
former  provides  the  better  fit,  and  it  is  suggested  that  it  be  used  if  temperatures  above  329  K  are 
not  eneountered. 
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The  strain-hardening  terms,  A,  B,  and  n,  ean  be  fit  to  the  stress-strain  eurves  of  figure  3,  under 
the  assumption  that  the  high-rate  eurves  are  adiabatie  and  the  low-rate  eurves  are  isothermal. 

This  is  done  in  an  ad  hoe  manner  to  roughly  mateh  the  data.  To  fit  the  adiabatie  eurves,  it  is 
assumed  that  ev,  the  speeifie  heat  eapaeity,  is  1160  J/kg-K,  based  on  an  estimate  used  by 
Raftenberg  et  al.  (7).  The  eonversion  of  plastie  work  to  heat,  |3,  is  assumed  to  be  0.9.  *  Finally, 
the  parameters  are  arbitrarily  adjusted  sueh  that  the  overall  fit  of  the  model  adequately  matches 
the  experimental  data  of  figures  3,  4,  6-8,  simultaneously,  using  the  actual  conditions  at  which 
those  graphs  present  the  data  (e.g.,  actual  plastic  strain  values,  adiabatic/isothermal  conditions). 
More  sophisticated  fitting  procedures  can  be  developed  but  are  hardly  justifiable  as  they  do  not 
offer  any  real  improvement  to  the  quality  of  the  fit  obtained. 

Four  sets  of  final  parameters  are  given  in  table  1,  for  each  model  over  the  two  temperature  ranges 
mentioned  above.  Comparisons  of  selected  model  fits  against  the  experimental  data  are  given  in 
figures  6,  8,  9,  and  10.  In  general,  the  difference  between  the  quality  of  the  fits  between  the  JC 
and  MJC  is  small,  with  the  MJC  providing  a  slightly  superior  fit  to  the  data  (see  figure  6)  in  a 
least-squares-residuals  sense.  Because  there  is  no  rate  dependence  of  the  strain-hardening  terms, 
neither  model  is  able  to  capture  the  strain-hardening  effect  at  both  high  (>600/s,  figure  9)  and 
low  (<0.1/s,  figure  10)  rates.  This  is  obvious  in  the  figures.  Because  it  is  anticipated  that  these 
models  will  be  applied  to  high-rate  scenarios,  the  data  is  fit  to  the  high-rate  data  essentially 
ignoring  the  strain  hardening  observed  at  the  low  rates.  Users  are  cautioned  against  using  these 
parameter  sets  if  their  problems  involve  significant  deformation  at  low  rates. 


Table  1.  JC  and  MJC  constants  for  Al/PTFE.  The  fits  are  performed  over  two  temperature  ranges. 


Parameter 

295  K  <  T  <  329  K 

295K<T<351K 

Unit  of 
Measurement 

JC 

MJC 

JC 

MJC 

A 

22.71 

22.71 

22.71 

22.71 

MPa 

B 

160 

160.1 

160 

160 

MPa 

n 

1.8 

1.8 

1.8 

1.8 

— 

c 

0.0339 

Not  used 

0.0339 

Not  used 

— 

k 

Not  used 

0.0324 

Not  used 

0.0324 

— 

(de/dt)o 

1 

1 

1 

1 

m 

1 

1 

0.707 

0.707 

— 

417 

417 

541 

541 

K 

Tr 

295 

295 

295 

295 

K 

P 

2290 

2290 

2290 

2290 

kg/in^ 

Cv 

1161 

1161 

1161 

1161 

J/kg-K 

Note:  Tm  as  used  here  is  simply  a  free-parameter  in  the  fitting  process  and  is  not  the  melting  temperature  of  the  material. 


Rae  and  Dattelbaum  (8)  measured  the  beta  factor  for  pure  PTFE  and  found  it  to  be  in  the  range  from  0.8  to  0.6  at  room 
temperature  and  at  a  strain  rate  of  10/s. 
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Figure  9.  MJC  model  predictions,  assuming  adiabatic  conditions,  compared  to  the  high-rate 
(SFIPB)  data.  The  JC  model  predictions  are  essentially  identical. 


true  strain 


Figure  10.  MJC  model  predictions,  assuming  isothermal  deformations,  compared  to  the  low-rate 
experiments.  The  JC  model  predictions  are  essentially  identical. 
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5.2  Zerilli-Armstrong  Model  for  Polymers 

Zerilli  and  Armstrong  have  developed  constitutive  models  for  metals  based  on  thermally 
activated  dislocation  motion  (P,  10).  Recently,  they  have  adapted  the  basic  framework  to  apply 
to  viscoplastic  deformation  of  polymers  {11, 12).  They  have  used  it  to  model  PTFE  (11),  and 
Cai  et  al.  (13)  have  used  it  to  model  a  composite  mixture  of  PTFE-Al-W.  Based  on  the  success 
of  these  applications,  it  was  decided  to  use  this  model  to  represent  the  present  data.  Although  the 
model  is  physically  based,  the  approach  taken  here  is  to  treat  it  as  a  curve  fit,  i.e.,  the  values 
presented  below  may  not  necessarily  be  realistic  beyond  the  observation  that  they  match  the 
existing  data. 


The  Zerilli-Armstrong  equation  for  polymers  relates  effective  stress,  cr,  to  effective  plastic  strain, 
£p,  temperature,  T,  and  pressure,  p,  according  to 


a  =  Be-P^  +  B 


1-e' 


,-aT 


(4) 


where 


and 


P  =  -  A  In  A’ 

(5) 

«  =  ctp  -  In  , 

(6) 

(0  =  0)^  Xne^  +0)pP, 

(V) 

b  =  b^(i+b,,p)“', 

(8) 

b.  =  b,,.('+b,„p)‘*'. 

(9) 

There  are  13  parameters  that  need  to  be  determined.  Due  to  the  relatively  limited  set  of  data  and 
the  complexity  of  the  model,  a  number  of  parameters  are  set  to  be  zero.  Foremost,  all  parameters 
leading  to  a  pressure  dependent  yield  stress  are  set  to  zero.  This  is  not  to  imply  that  the  material 
has  no  pressure  dependence,  the  data  is  simply  not  available.*  Thus,  Bpb  =  cOp  =  Bopn  =  Bpn  =  Bopb 
=  0,  and  the  model  simplifies  to 


ct  =  B  +  R 

pa  Opa-' 


h  _  I  ,  .  1 

l-I-  ^ 

[co^+co,\ne^) 


(10) 


Assuming  £p  =  0  and  T  =  295  K,  the  right-hand  term  vanishes  and  a  power-law  least-squares  fit 
to  the  rate  data  (figure  6)  can  be  made. 


cr  = 


(7  =  ke^ 


Ar 


(11) 


PTFE  has  been  shown  to  have  a  pressure  dependent  flow  stress  (.10). 
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The  constant  k  and  |3i  are  determined  by  the  curve  fit.  Through  k,  a  relationship  is  established 
between  |3o  and  Bpa. 


Po  = 


f 


V 


(12) 


Using  this  relationship  and  assuming  £  =  4000/s  and  £p  =  0,  Bpa  can  be  determined  by  that  value 
which  gives  an  optimum  fit  to  the  temperature  data  of  figure  8.  Equation  12  then  gives  |3o,  and 
the  first  term  on  the  right-hand  side  of  equation  10  is  determined. 


The  remaining  term  deals  with  strain  hardening.  Aside  from  the  behavior  seen  in  figure  7  at 
351  K,  which  may  be  due  to  a  phase  change,  there  is  no  convincing  data  that  the  strain  hardening 
changes  with  temperature,  so  set  ao  =  ai  =  0.*  This  leaves  three  terms  to  be  determined;  (Oa,  (Ob, 
and  Bopa.  A  least  absolute  value  best  fif*^  is  made  to  the  stress-strain  data  shown  in  figure  3, 
assuming  the  high-rate  curves  are  isothermal,  i.e.,  the  data  are  fit  to 


cr  =  5 

pa  0  pa ' 


1-e 


ln£p  )fp 


(ty,+ty,lnfj  ’ 


(13) 


to  obtain  estimates  for  cOa,  (Ob,  and  Bopa.  This  is  technically  incorrect,  but  fortunately  the 
differences  between  adiabatic  and  isothermal  behaviors  for  these  experiments  are  not  drastically 
different,  and  thus  the  fitting  process  is  simplified. 

As  with  the  Johnson-Cook  fits,  the  resulting  parameter  set  is  adjusted  in  an  ad  hoc  manner  to 
obtain  the  “best”  simultaneous  fits  to  the  isothermal  and  adiabatic  stress-strain  curves,  the  rate 
hardening,  and  the  thermal  softening.  The  resulting  optimum  parameters  are  listed  in  table  2, 
and  the  quality  of  the  fit  is  shown  against  the  experimental  data  in  figures  6,  8,  11,  and  12.  To 
facilitate  comparisons,  the  parameter  sets  of  Cai  et  al.  (75)  for  PTFE-Al-W,  and  Zerilli  and 
Armstrong  (77)  for  PTEE  and  also  listed  in  the  table. 


6.  Conclusion 


The  experiments  conducted  in  this  program  have  been  used  to  populate  the  JC,  MJC,  and  ZA  for 
polymers  constitutive  models.  These  models  address  the  inert  behavior  of  the  material  only.  It  is 
worth  noting  that  in  none  of  these  experiments  was  any  evidence  of  a  reaction,  either  during  the 
experiment  or  post-mortem,  observed. 


* 

It  was  hoped  that  adjustment  of  Uq  and  tti  would  help  produce  a  better  match  between  the  model  and  the  temperature  tests  at 
larger  deformations.  However,  this  did  not  appear  to  be  the  case  so  these  two  parameters  were  left  at  zero  values, 
numerical  technique  is  used. 
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Table  2.  The  ZA  constants  for  Al/PTFE.  The  values  detemined  by  Cai  et  al.  (13)  and  Zerilli  and 
Armstrong  (11)  for  PTFE-Al-W  and  PTFE  are  included  for  comparison. 


Parameter 

RM-4 

PTFE-Al-W 

ZA  PTFE 

Unit  of 
Measurement 

Po 

0.011672 

0.020100 

0.020100 

1/K 

Pi 

0.000139 

0.000264 

0.000264 

1/K 

ao 

0.000000 

0.004780 

0.004780 

1/K 

ai 

0.000000 

0.000050 

0.000050 

1/K 

COa 

-3.000 

-2000 

-3.600 

— 

(flb 

-0.500 

-0.625 

-0.625 

— 

Md 

0.000 

-0.031 

-0.040 

1/MPa 

Bpa 

550 

4016 

4016 

MPa 

Bpb 

0.000 

0.020 

0.020 

1/MPa 

Bpn 

0.000 

0.714 

0.714 

— 

^Ooa 

25.0 

72.4 

72.4 

MPa 

^Oob 

0.000 

0.022 

0.022 

1/MPa 

^Onn 

0.000 

0.500 

0.500 

— 

Figure  11.  Comparison  of  the  ZA  model  to  the  SFIPB  experiments  assuming  adiabatic  conditions. 
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true  strain 


Figure  12.  Comparison  of  the  ZA  model  to  the  low-rate  experiments  assuming  isothermal  conditions. 

Note  the  great  improvement  over  the  JC  fits  to  the  same  data,  figure  10. 

It  is  believed  that  these  models  are  suitable  for  simulations  of  low-speed  impacts  (e.g.,  the 
cylinder  impacts  of  Casern  and  Raftenberg  [74]  and  Mock  and  Holt  [75]).  However,  readers 
shown  remain  aware  of  the  comments  and  caveats  made  above,  especially  the  inability  of  the  JC 
models  to  capture  the  strain  hardening  effects  over  the  full  range  of  strain-rates  investigated.  It  is 
suggested  that  the  JC  model  could  be  made  to  capture  this  effect  by  making  the  strain  hardening 
parameters  suitable  functions  of  strain-rate.  However,  it  is  recommended  that  the  ZA  model, 
which  is  more  firmly  rooted  by  a  physical  foundation,  be  used  in  place  of  the  JC  fits  whenever 
possible.  As  with  any  model,  readers  should  exercise  caution  when  using  these  models  outside 
of  their  calibrated  ranges.  Note  that  the  elastic  behavior  of  the  material  has  not  been  characterized. 
Any  reader  who  implements  this  model  will  need  to  address  this.  One  should  also  take  care  when 
using  these  models  to  represent  Al-PTFE  mixtures  that  are  not  manufactured  using  this  same 
process.  For  example,  there  is  a  considerable  difference  between  this  material  behavior  and  that 
measured  in  Raftenberg  (7),  although  superficially  both  are  basically  the  same  mixture. 

Finally,  failure  is  not  addressed  in  this  work.  When  failure  was  observed  in  these  tests,  it 
appeared  to  be  primarily  due  to  ductile  fracture.  It  is  not  believed  that  this  material  exhibits 
shear  banding,  and  the  mechanical  properties  listed  in  this  work  support  this  (in  terms  of  the 
relative  effects  of  strain  and  strain-rate  hardening  and  thermal  softening).  Some  fracture  work 
has  been  preformed  on  Al/PTFF  composites;  refer  to  Joyce  and  Joyce  (76). 
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Appendix.  Strain  Measures  Defined 


This  appendix  clarifies  the  definitions  of  engineering  stress  and  strain  and  true  stress  and  strain 
as  used  in  this  report.  The  unqualified  use  of  any  of  these  terms  implies  positive  values  in 
compression. 


Engineering  stress,  s,  is  defined  as 


(A-1) 


where  P  is  the  compressive  load  applied  to  the  specimen  and  Ao  is  its  original  cross-sectional 


True  stress,  a,  is  defined  as 


P 


(A-2) 


where  A  is  the  instantaneous  cross-sectional  area  of  the  specimen,  i.e.,  the  area  at  any  instant  in 
time  as  the  specimen  deforms. 


Engineering  strain,  e,  is  defined  by 


(A-3) 


where  5  is  the  amount  the  specimen  is  compressed  and  Eq  is  the  initial  length  of  the  specimen. 
Note  that  E,  the  instantaneous  length,  can  then  be  found  by  E  =  Eg  -  5 . 


True  strain,  £,  is  defined  by 


E 


(A-4) 


Here,  dL  is  an  increment  in  specimen  length,  i.e.,  the  negative  sign  is  needed  to  make  this 
positive  in  compression.  This  can  be  integrated  for  the  conditions  of  the  compression  test  to  give 


£  =  -ln  —  . 

E  A  , 


(A-5) 


By  equation  (A-3),  true  strain  can  be  related  to  engineering  strain  by 


£  =  -  ln(l  -  e)  . 


(A-6) 
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If  the  deformation  of  a  specimen  is  incompressible,  and  if  the  initially  cylindrical  specimen 
remains  cylindrical  during  the  deformation  process,*  the  instantaneous  volume  is  equal  to  the 
initial  volume  at  any  instant,  i.e., 

d^L  =  do"Lo,  (A-7) 


where  do  and  d  are  the  initial  and  instantaneous  diameters  of  the  specimen.  Using  this 
relationship,  true  stress  can  be  determined  from  the  engineering  quantities  by 

a  =  s(l-e).  (A-8) 


Furthermore,  the  tensile  radial  engineering  strain  of  the  specimen,  defined  as 


d-d„ 


e,  = 


(A-9) 


can  be  related  to  the  axial  engineering  strain  of  the  specimen  by 

(A-10) 

Vl  -  e 

This  applies  only  to  incompressible  behavior.  It  is  this  equation  that  was  plotted  in  figure  2  of 
this  report. 


This  is  only  approximately  true  in  compression  testing  because  of  “barreling”  that  occurs  due  to  friction  at  the  platens. 
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